Motion of charged particles near weakly magnetized Schwarzschild black hole 
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We study motion of a charged particle in the vicinity of a weakly magnetized Schwarzschild black 
hole and focus on its bounded trajectories lying in the black hole equatorial plane. If the Lorentz 
force, acting on the particle, is directed outward from the black hole, there exist two qualitatively 
different types of trajectories, one is a curly motion and another one is a trajectory without curls. 
We calculated the critical value of the magnetic field for the transition between these two types. If 
the magnetic field is greater than the critical one, for fixed values of the particle energy and angular 
momentum, the bounded trajectory has curls. The curls appear as a result of the gravitational drift. 
The greater the value of the magnetic field, the larger is the number of curls. We constructed an 
approximate analytical solution for a bounded trajectory and found the gravitational drift velocity 
of its guiding center. 



PACS numbers: 04.70.Bw, 04.70.-s, 04.25.-g 



I. INTRODUCTION 

There exist both theoretical and experimental indica- 
tions that a magnetic field must be present in the vicinity 
of black holes. A regular magnetic field can exist near 
a black hole surrounded by conducting matter (plasma), 
e.g., if the black hole has an accretion disk. The magnetic 
field near a stellar mass black hole may contain a contri- 
bution from the original magnetic field of the collapsed 
progenitor star. The dynamo mechanism in the plasma 
of the accretion disk might generate a regular magnetic 
field inside the disk. Such a field cannot "cross" the con- 
ducting plasma region and is trapped in the vicinity of 
the black hole (see, e.g., discussion in [l|). 

Stellar mass and supermassive black holes often have 
jets, that is the collimated fluxes of relativistic plasma, 
ft is generally believed that the MHD of the plasma in 
strong magnetic and gravitational fields of the black holes 
would allow one to understand formation and energetics 
of the black holes jets [l[ . Magnetic field in the vicinity of 
the black holes might play an important role in transfer 
the energy from the accretion disc to jets. Existence of a 
regular magnetic field near the black holes is also required 
for the proper collimation of the plasma in the jets. 

The simplest mechanism for extraction of the rotation 
energy from black holes in the presence of magnetic field 
was proposed by Blandford and Znajek m. This mecha- 
nism is discussed in detail in the nice book [3[ . Interesting 
results of the numerical simulations, demonstrating jet 
formation by the plasma in strong magnetic and gravita- 
tional fields of the black holes, are presented in [J] (see 
also references therein). There arc some observational 
evidences of the existence of magnetic field around back 
holes and its effect on the dynamics of their accretion 
disks (see, e.g., [!, Hi)- 
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In this paper, we study motion of charged particles in 
magnetized black holes. To make our presentation more 
concrete, we shall use the estimations for the magnetic 
field given in Q • Namely, the characteristic scales of the 
magnetic field B arc of the order of B\ ~ 10 8 G near 
horizon of a stellar mass (Mi ~ 10M Q ) black hole, and 
of the order of B2 ~ 10 4 G near horizon of a supermas- 
sive (M2 ~ 10 9 M Q ) black hole. The value of B2 is also 
discussed in connection with the Blandford-Znajek mech- 
anism (see, e.g., [H, H|). It should be emphasized that the 
aim of the present paper is not to discuss special astro- 
physical objects, but to study interesting features of the 
dynamics of a charged particle motion in the presence 
of both the gravitational field of such black holes and 
the magnetic field in their vicinity. We use the estimates 
above to describe only the domain of the physical pa- 
rameters, which characterize our system, and to define 
the validity of the adopted approximations. 

Let us notice that both the fields B\ and B2 are weak 
in the following sense: The space-time local curvature 
created by the magnetic field B is of the order of GB 2 /c 4 . 
It is comparable to the space-time curvature near a black 
hole of mass M only if 
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For a black hole of mass M this condition holds if 
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Evidently, the quantities B\.i for both the stellar mass 
and supermassive black holes presented above are much 
smaller than the corresponding Bm- This means that 
for our problem the field B can be considered as a test 
field in the given gravitational background. Such a mag- 
netic field practically does not affect motion of neutral 
particles. 

However, for charged particles, the acceleration in- 
duced by the Lorentz force can be large. The acceler- 
ation is of the order of qB/(mc). Thus, the 'weakness' of 
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the magnetic field B is compensated by the large value 
of the charge-to- mass ratio q/m, which, for example for 
electrons, is e/m e w 5.2728 x 10 17 g _1/2 cm 3 / 2 s- 1 . 

In this paper, we consider a Schwarzschild black hole 
immersed into magnetic field, which is homogeneous at 
the spatial infinity, where it is directed along the 'verti- 
cal' z-axis. Close to the black hole this field is strongly 
affected by its gravity. Such a model is a good approx- 
imation for more realistic magnetic fields generated by 
currents in conducting accretion disk, provided the size 
of the black hole is much smaller than the size of the 
disk (see, e.g., discussion in Q). In our analysis, we shall 
neglect mutual interaction of charged particles moving 
around the black hole, i.e., we consider the approxima- 
tion of a diluted accretion disk. We also restrict our- 
selves to the motion in the equatorial plane of the black 
hole, which is orthogonal to the direction of the magnetic 
field. Similar approximations were considered in the pa- 
pers 0,113. 

In the presence of the external magnetic field the 
Lorentz force, acting on a charged particle, depends on 
the direction of the particle motion. Namely, the pres- 
ence of the magnetic field breaks the discrete reflection 
symmetry in the azimuthal direction, so that motion in 
clockwise and counter-clockwise directions is not equiva- 
lent. Many important results on a charged particle mo- 
tion near a magnetized black hole can be found in the 
interesting review by Aliev and Gal'tsov [lj| ■ A study 
of the marginally stable circular orbits around a magne- 
tized Kerr black hole can be found in ll(. Motion of 
charged particles in the dipolar magnetic field around a 
static black hole was studied in [l^Hl4j . and around a 
rotating black hole in (l5| . 

To compare the relative strength of the magnetic and 
gravitational forces acting on a charged particle moving 
around the black hole let us make the following simple 
estimations: In a flat space-time (in the absence of grav- 
ity) a particle with charge q and the rest mass mass m 
in the magnetic field B has the characteristic cyclotron 
frequency 



a..= 



\qB\c 
E ' 



(3) 



where E = 7 mc 2 is energy of the particle. Let us com- 
pare this frequency with the Keplerian frequency of a 
particle orbiting a black hole of mass M, 



K 



3/2^2 



(4) 



Here r g = 2MG/c 2 is the black hole gravitational radius. 
For r ~ r g and 7 *~ 1 the ratio of these frequencies Q, c /Q,k 
is of the order of 



b = 



qBMG 



(5) 



One can use this dimensionlcss quantity to characterize 
the relative strength of the magnetic and gravitational 



forces acting on a charged particle moving near the black 
hole. For the motion far from the black hole, the ratio 
f2 c /f2if contains the additional factor (r/r g ) 3 ^ 2 . 

For a proton of the mass m p and the charge e moving 
near a magnetized stellar mass or supcrmassivc black hole 
wc have 



eB\GM\ 



4.7180 x 10 7 



B 2 M 2 
B-lMi 



b\ = 10% 



(0) 
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respectively. Both the quantities are large. For electrons 
their value is larger by the factor m p /m e ~ 1836. This 
indicates that the effect of the magnetic field on a charged 
particle motion is not weak. In a general case, such a 
magnetic field essentially modifies the motion of charged 
particles. For example, due to this effect the innermost 
stable circular orbits around a Schwarzschild black hole 
for proper orientation of the particle rotation are shifted 
towards its horizon (see, e.g., [9I4TTI. [l6| ) . 

The aim of this paper is to study the dynamics of a 
charged particle moving near magnetized black holes. In 
particular, we shall demonstrate that in the presence of 
the magnetic field, there exist two qualitatively different 
types of the bounded orbits. In the first type, the Keple- 
rian motion is modulated by the cyclotron revolution, so 
that the particle trajectory has a curly-type structure at 
small scales. In the second type the particle trajectory 
has no curls. We study the transition between these two 
types and calculate the critical value of the corresponding 
magnetic field. A combined action of gravity and elec- 
tromagnetism shifts the Keplerian radius and generates 
the drift of the guiding center of the particle trajectory. 

This paper is organized as follows: In Sec. II we 
present equations of motion of a charged particle moving 
around a weakly magnetized Schwarzschild black hole. 
In Sec. Ill we study motion of a charged particle in a 
uniform magnetic and weak gravitational fields. Section 
IV contains an analysis of the effective potential corre- 
sponding to a charged particle moving in the equatorial 
plane of the black hole. In Sec. V we study the innermost 
stable circular orbits. Section VI contains an analysis of 
bounded trajectories. In Sec. VII we present an approx- 
imate solution for bounded motion of a charged particle. 
We summarize our results in Section VIII. In this paper 
wc use the sign conventions adopted in [It} and units 
where c = 1 . 



II. EQUATIONS OF MOTION 

A. Magnetized black hole 

We shall study a charged particle motion in the vicinity 
of a Schwarzschild black hole of mass M in the presence 
of an axisymmetric and uniform at the spatial infinity 
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magnetic field B. The Schwarzschild metric reads 



ds 2 = -fdt 2 + f^dr 2 



2 duj 2 



f = 1 - ^ , (8) 
r 



where r g = 2GM and cLu 2 = d6 2 + sin 2 6d(b 2 . The com- 
muting Killing vectors £™ = d/dt and ^) = 9/90 gen- 
erate time translations and rotations around the symme- 
try axis, respectively. 

Because the metric ((SJ) is Ricci flat, the Killing vectors 
obey the equation 



0. 



(9) 



This equation coincides with the Maxwell equation for 
a 4-potential A a in the Lorentz gauge A a . a = 0. The 
special choice 



-r 

2 H<P) 



(10) 



corresponds to a test magnetic field, which is homoge- 
neous at the spatial infinity where it has the strength 
B (see, e.g., [HI EH). l n what follows, we assume that 
B > 0. The electric 4-potential (fTU|) is invariant with re- 
spect to the isometries corresponding to the Killing vec- 
tors, i.e., 



{L^A) a = A a ^» + A b i h a = . 
The magnetic field measured in the rest frame is 



B a 



2 '~abcd r 



\£;-\ 1/2 



£0123 — 
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(ii) 
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Here F ab = A bM — A a , b is the electromagnetic field tensor. 
The magnetic field corresponding to the 4-potcntial (fTU)) 
is 



r \ 1/2 
B a d a = B (1— — 



cos ( 



d_ 

a7 



d 



r 86 



(13) 



In the cylindrical coordinates (p, z,(j)) at the spatial in- 
finity it is directed along the z-axis. The magnetic field 
(|13[) coincides with the dipole approximation of the mag- 
netic field of a current loop of the radius r = a around 
a Schwarzschild black hole One can show that this 
approximation is valid in the region r g < r <C a. 



B. Dynamical equations 



Dynamical equation for a charged particle motion is 

mu a = qF a b u b , (14) 

where u a is the particle 4- velocity, u a u a = — 1, q and m 
are its charge and mass, respectively. Here and in what 
follows, we denote by the over dot the derivative d/dr 
with respect to the proper time r. For the motion in 
the magnetized black hole (© and ([TO")) ) there exist two 



conserved quantities associated with the Killing vectors: 
the energy E > and the generalized angular momentum 

L e (—00, +oo), 



L 



br 2 sin 2 I 



qB 
2 



9 ' 2 

r sin i 



(15) 
(16) 



Here P a = m u a + qA a is the generalized 4-momentum of 
the particle. 

It is easy to check that the ^-component of Eq. fTJJ 
allows for a solution 9 = n/2. This is a motion in the 
equatorial plane of the black hole, which is orthogonal 
to the magnetic field. We restrict ourselves to this type 
of motion, for which the conserved quantities (|15p and 
(ITBl arc sufficient for the complete integrability of the 
dynamical equations. 

Let us denote 



m 



m 



2m ' 



(17) 



where £ and C are the specific energy and the specific 
angular momentum. The r-component of the dynamical 
equation (fT4")) for the planar orbits is 



B 



ilTg 

2r 



B 1 



JL 
2r 



'9 

2r 2 



(18) 

The first integral of this equation, which follows from 
Eqs. ([15]). (fT6|) . and dTSj) is given by 



= £-(!- — 



1 



(19) 



The equations for the azimuthal angle <f> and the time t 
are 



= — -B, t = £\\-S- 



(20) 



Equations (fT8|) -([20 |) are invariant under the following 
transformations: 



(21) 



Thus, without loss of the generality, one can assume that 
the charge q (and hence B) is positive. For a particle 
with a negative charge it is sufficient to make the trans- 
formation (|2"Tj) . 

The geodesic motion in the Schwarzschild space-time 
UHJ is well studied (see, e.g., [l3, [T3, Here we study 
how this motion is modified in the presence of the mag- 
netic field (|T3"|) . Our problem is reduced to the set of three 
first-order differential equations (fT9|) and ([20]) containing 
two independent parameters £ and £, which are integrals 
of motion. The solution to the problem depends on the 
parameters: ro, </>o, arid to. They define the initial po- 
sition of the charged particle. The radial equation (fT9|) 
can be integrated independently. Substituting its solu- 
tion into Eq. (|20l) one can obtain the angular variable <fi 
and the time t by integration. 
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III. WEAK GRAVITATIONAL FIELD 
A. Flat space-time limit 

Before we study the motion of a particle near the black 
hole, it is instructive to consider an approximation of a 
weak gravitational field. In the absence of gravity, when 
r g = 0, equations (|T5|) and (f^Uj) arc greatly simplified and 
take the following form: 

r 2 =£ 2 -l-r 2 (^§-s) 2 , (22) 

4>=^-B, i = £. (23) 

The factor £ in a flat space-time coincides with the 
Lorentz 7-factor. The solution to equations (|22[) and (|23|) 




(a) (b) 

FIG. 1. Motion of a charged particle in a uniform magnetic 
field in a flat space-time. Points A and B define the maximal 
and the minimal distances of the particle from the coordinate 
origin O. (a): C < 0. (b): £ > 0. Points C and D are the 
turning points where = 0. 

is well known. The particle moves along a circle of the 
radius 

r c = ^\£\/B, (24) 

where C < is defined with respect to the center of the 
circle O' . The corresponding cyclotron frequency ([3]) and 
the Lorentz 7-factor are 

2B 

O c = — , 8= Vl + 4B|£|. (25) 

The general solution to equations (f2"2"j) and (|23p. which 
contains two additional constants r$ and </>0j can be easily 
obtained by coordinate transformation sending the center 
of the circle to an arbitrary point O on the plane z = 0. 
Figure [TJa) illustrates the case when the coordinate ori- 
gin is located inside the circle. More interesting for us 
case is shown in Fig.QJb), where the coordinate origin is 
located outside the circle. For such choice of the coordi- 
nates the motion in the r-direction is bounded, so that 
Tmin < r < ''max, and the angle <j> changes between the 
two values 0_ and <p + . Motions in the r- and (^-direction 
are correlated and have the same period, so that the re- 
sulting trajectory is a circle. 



B. Approximation of a weak gravitational field 

In the presence of the black hole there is a preferable 
choice of the origin O of a coordinate system, namely the 
black hole location. Thus, the symmetry of the flat space- 
time solution is broken. A black hole also is a source of 
the additional (gravitational) force acting on the parti- 
cle. Let us now discuss the motion of the particle in the 
homogeneous magnetic field, in the approximation when 
this force is weak. 

The corresponding equations can be obtained from 
equations (JTSJ and (|20[) . assuming that r g /r -C 1. In 
the leading order we have t « £t, and 

f=^-B 2 r- 9l (26) 

^^I- B - ( 27 ) 

Here g = r g /(2r 2 ) is the Newtonian gravitational force. 
This force is directed along the radius toward the black 
hole, and is orthogonal to the magnetic field. 

Take a point (ro, 4>o), such that r g <C Tq, and introduce 
the local Cartesian coordinates (x, y) near it, so that y 
is directed along the radius vector at the point, and x is 
directed "clockwise". For x -C and y ^ r , which cor- 
responds to short time intervals At and strong magnetic 
field B, we have 

r « r + y , cf> 4>o - x/r . (28) 

Since g contains the small factor r g /ro, it is enough to 
keep the leading term 

So = X ■ (29) 
z 'o 

We substitute expansion (f2"5| into equations (|2"o) and 
(|27|l. and keep the terms linear in x and y. 
In the zero order, these equations give 

C = Br 2 . (30) 

The first order terms give the following dynamical equa- 
tions: 

y = -n 2 y-g , (31) 
x = Qy . (32) 

Here f2 = 2B = Vl c £ is the relativistic cyclotron fre- 
quency. 

The solution to equations [) and (|32p corresponding 
to the initial conditions x(Q) = and y(0) = a — go/fl 2 
is 

2/(r) = aco S (Or)-^, (33) 

x(t) = asin(fh-) - (go/tyr . (34) 

The line on the (x, y) plane described by this solution is 
called a trochoid. For go = this solution describes the 
motion along a circle with the cyclotron frequency fl c . 
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FIG. 2. Motion of a charged particle in the mutually orthog- 
onal uniform magnetic and weak gravitational fields. Dashed 
line illustrates the motion of the guiding center. 



In the presence of gravity the center of the circle moves 
in the negative ^-direction with the velocity 



v = r ( ( />) = -. 



(35) 



The velocity V defined with respect to the rest frame 
differs by the Lorentz factor t = £ , 



V=-= 90 



£ n c £ 2 



(36) 



IV. CHARGED PARTICLE IN THE 
SCHWARZSCHILD SPACE-TIME 

A. Dynamical equations 

We return to our main problem, motion of a charged 
particle near a magnetized black hole. It is convenient 
introduce the following dimensionless quantities: 



T = 



P = 



t = 



C 



(39) 



and write B = b/r g (see Eqs. © and (fT7]l). Then, the dy- 
namical equations (fTi?|) and (|20l) take the following form: 



Here 



da 



U = 



dp 
da 
(_ 

^2 



1-i 

P 



dT 

da 



£p 



1 



{t-bp 2 f 



(40) 
(41) 

(42) 



is the effective potential. According to the adopted con- 
vention, we have b > 0. 

The parameter £ can be positive or negative, £ = ±|£|. 
For t > (sign +) the Lorentz force, acting on a charged 
particle, is repulsive, i.e., it is directed outward from the 
black hole, and for t < (sign — ) the Lorentz force is 
attractive, i.e., it is directed toward the black holc0. 



For a = ga/fl 2 the solution takes the form 

J/(T) = §[cos(Or)-l], (37) 
x(t) = ^[sin(Or)-Or]. (38) 

The line on the (x, y) plane described by this solution 
is called a cycloid. It separates two different types of 
motion, with curls (for a > go /ft 2 ) and without (for 
a < go /ft 2 ) (see Figure [2|). Such a curly- type structure is 
characteristic for a charged particle motion in the pres- 
ence of both magnetic and non- magnetic forces (see, e.g., 
IU, Chapter 2). 

The obtained solution has a simple interpretation. In a 
frame moving with the velocity (|51)|) the magnetic field B 
'generates' the electric field E y = ~fVB, which is orthog- 
onal to both the field B and the velocity V. This elec- 
tric field, directed along the y-axis, acts on the charge 
q with the force qE. The velocity V is defined by the 
condition, that this force exactly compensates the grav- 
itational force mgo- Thus, the motion of a charged par- 
ticle due to the gravitational force exerted on it in the 
direction orthogonal to the magnetic field is analogous to 
the motion of the particle in the constant and mutually 
orthogonal electric and magnetic fields (see, e.g., [22j]). 



B. Effective potential 

Let us study properties of the effective potential. We 
assume that the mass of the black hole and the strength 
of the magnetic field are fixed. Thus, for a given charged 
particle the parameter b is also fixed, and the effective 
potential is a function of two variables, I and p. We 
focus on the motion in the black hole exterior, where 
p > 1. Expression (|4"2"j) shows that the effective potential 
is positive in that region. It vanishes at the black hole 
horizon, p = 1, and grows as b 2 p 2 for p — > +oo. The 
latter property implies that the particle never reaches 
the spatial infinity, i.e., its motion is always finite. 

To study the characteristic features of U we use its first 
and second derivatives 

U. p = y{2b 2 p b - 6V ~ 2£bp 2 + p 2 - 2£ 2 p + 3l 2 ) , (43) 
U, pp = 2- (b 2 p 5 -p 2 + 2£b p 2 + U 2 p - U 2 ) . (44) 

1 According to the terminology used in the papers 0] and lldl . 
"Larmor" rotation corresponds to the Lorentz force pointing to- 
ward the black hole, while "anti-Larmor" rotation corresponds 
to the Lorentz force pointing in the opposite direction. 
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Extrema of the effective potential are defined by the 
equation U iP = 0. The extremal points are the roots of 
the fifth order in p polynomial in the brackets of expres- 
sion (|43|) . Such a polynomial may have five real valued 
roots. We shall show now that because of the specific 
structure of this polynomial, it has at most two real val- 
ued roots in the region p > 1. 

For this purpose we write the equation U. p — in the 
following equivalent form: P(p) = Q(p), where 

2 



4- 



P(p)=b 2 p%2p-1) 



Q(p) = 2£bp 2 



2£ 2 p 



P 
- 3£ 2 



(45) 
(46) 



2- 



To find a solution to this equation corresponding to the 
black hole exterior, it is sufficient to find points of the 
intersection of the curves y = P{p) and y = Q(p) in the 
interval p 6 [l,oo). Note that in this interval P(p) is 
a positive and monotonically growing function, without 
points of convolution. The second function represents a 
parabola, which can either intersect the curve y = P{p) 
at two points, or do not intersect it at all. Note that the 
case when there is only one intersection in the black hole 
exterior is excluded for the following reason: In such a 
case the potential would have only one extremum (max- 
imum or minimum), what is impossible for a positive 
function vanishing at p = and growing at infinity. In 
the limiting case these curves touch each other at one 
point, where P p = Q p . It is easy to check that at this 



1- 



u 
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£ 2 

m ax. |. 
f ' 
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/ f *\ 
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pi 
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1 Pn 
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4 



point U , pp 



= 0. 



In the former case the intersection of the curves occurs 
at two points, p = p max and p = p min > p max . The 
maximum of U is at p = /9 max , while its minimum is at 

P = Pr, 



FIG. 3. Effective potential for I > and b = 1/2. Curve 1 
corresponds to £ sd 1.18. There are no circular orbits. Curve 
2 corresponds to t sd 2.07. There is the innermost stable cir- 
cular orbit defined by p = p+ w 1.59. Curve 3 corresponds to 
£ ss 3.22. There are unstable and stable circular orbits corre- 
sponding to the minimum and the maximum of the effective 
potential. For the bounded trajectory corresponding to the 
energy S+ one has p € [pi,p 2 ]. 



cyclotron and Keplerian oscillations due to the gravi- 
tational and Lorentz forces acting on a charged parti- 
cle. In the Schwarzschild space-time circular orbits cor- 



This means that at p = p min there is a stable rcsporidin g to rmin > 3rg are a i ways stab i ej i.e., ft2 > 0. 



circular orbit, while at p — p max there is an unstable 
one. The condition U tPP — singles out the innermost 
(or marginally) stable circular orbit. Figure [3] illustrates 
the behavior of the effective potential. 

Let us define the dimensionless parameter 

o I 



Marginally stable circular orbits correspond to f2 = C0. 

Let us mention another property of the effective po- 
tential, which we shall use later. We denote 



-U,pp{pmm) > , 



(47) 



which is a measure of the curvature of the effective po- 
tential at its minimum. It is equal to the square of the 
frequency of small oscillations about a circular orbit in 



p* = \fljb, 
then simple calculations give 

U p (p*) = b/£. 



(50) 



(51) 



the radial direction. One can use the equation U , 







This means that for bounded orbits and for positive £ 
one has 



to express £b in terms of other quantities. Substituting 
this expression into Eqs. (j42|) and (j4T|) . and eliminating 
£ one obtains the following expression for tu : 



P* > Pmin > Pn 



(52) 



*-'min(Pmin ^) 
^Pmin (Pmin ~ I) 2 



46 2 



"(Pn 



1) 



(48) 



where £ min = U(p m in)- This frequency can be written in 



the dimensional form as follows: 
3r 



n 2 K 



i 



i 



V. THE INNERMOST STABLE CIRCULAR 
ORBITS 

The equations U tP = and U iPP = 0, which determine 
the innermost stable circular orbits, can be presented as 



(49) 



where the Keplerian frequency Q,x is defined at r = r min . 
Such a relation was derived in @, E3 (for a more gen- 
eral case see [HI). It illustrates a combination of the 



2 The regions of existence and stability of the circular orbits with 
respect to small oscillations in the radial and vertical (9-) direc- 



tions where studied in 0. 
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(a) (b) 

FIG. 4. Dependance of the innermost stable circular orbits p = p± on b. For each value of b one has p+ < p_. For i> — > +00 
one has p+ — >• 1 and p_ — > (5 + \/l3)/4 w 2.15. 




0.0 0.1 0.2 0.3 0.4 0.0 0.05 0.1 0.15 0.2 0.25 

(a) (b) 

FIG. 5. Angular momenta t± for the innermost stable circular orbits as functions of 6. (a): £ > 0, the Lorentz force is repulsive, 
(b): £ < 0, the Lorentz force is attractive; the curve I = £-(b) is monotonically decreasing. For b = we have £± — ±\/3, what 
corresponds to a Schwarzschild black hole. For b — > +00 we have £± — > ±00. 



follows: 



&V(2p - 1) + P ~ 2£bp A - 2tp + 3r = , (53) 



0. 



(54) 



For a given 6 these equations allow one to find the pa- 
rameters p and I for such orbits. As usual for this kind of 
problems, it is much easier to obtain a solution in an im- 
plicit (parametric) form. An addition of equations (|53[) 
and (|54[) let us exclude the terms linear in I. Solving the 
obtained equation we derive 



±b 



4(3p±-l) 
(3 -Pi) 1 / 2 



1/2 



(55) 



The solution exists in the interval p± £ (1/3,3], where 
p± = 3 corresponds to b = 0. Since we study the black 
hole exterior, we have p± £ (1,3]. A substitution of 



Eq. (|55|) into either of equations (|53|) or (|54)l gives the 
following equation for p±: 



4/4-9p±+3± V(3p± - 1)(3 - p±)- 



3 - p-t 



= 0. (56) 



26 2 p| 

This equation allows one to express b in terms of p± 



b = 



2 P± Upl - 9p± + 3 ± y/(3 P± -l)(3-p±) 



1/2 ■ 



(57) 



Substituting this expression into Eq. (|55[) we obtain l± 
as a function of p±, 



±- 



^(Spi-l) 1 / 2 



v^2 (44 - 9p± + 3 ± V(3p±-l)(3-/9±; 



1/2 



(58) 
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The condition that b is real does not impose any new 
restrictions on p + , so that as above, 1 < p + < 3. For p_ 
it gives an additional restriction, which is (5 + vT3) /4 < 
p_ < 3. The position of the innermost stable circular 
orbit around a Schwarzschild black hole corresponds to 
b = 0, p± = 3, and i± = ±y/3. 

Figure @] shows plots of p± as functions of b. These 
plots demonstrate that if the magnetic field increases, 
the corresponding values of p + and p_ decrease. For a 
strong magnetic field and the repulsive Lorentz force, the 
radius of the innermost stable orbit may be very close to 
the gravitational radius [IJ (HI ■ Similar plots for t± as 
functions of b are shown in Figure [5] Let us mention 
that the curve I = £+(b) has the minimum i+ = (3 + 

y5) 2 /(8^2 + V5) w 1.67 at 6 = (V5 - 2) 3 / 2 /2 w 0.057. 
In the limit of the strong magnetic field one has 



da p 2 



(65) 



P+I&>>1 = l + w% +(D{b2) 



(59) 



4 36V136 2 



5 + ^13 41-11V13 ^ /t _ 4 , 

P-U»l = 7 + -„ /T77.0 +0 ( 6 )' ( 60 ) 



(61) 
(62) 



The leading terms in expressions (|59")l and (j6T))) were 
derived in |ll[ Il3j |. For stable circular orbits one has 
Pmin > P+ and I > £+, ii the Lorentz force is repulsive, 
and p m i n > p_ and I < ^_ , if the Lorentz force is attrac- 
tive. 



If the Lorentz force is attractive (£ < 0) the right hand 
side of this equation is always negative. The corre- 
sponding motion of a charged particle is in the clock- 
wise direction. This motion is modulated by oscillations 
in the radial direction. The corresponding trajectories 
have no curls. They are similar to bounded trajectories 
of test particles moving near a Schwarzschild black hole 
(see, e.g., [ill). The presence of the magnetic field just 
smoothly deforms them. 

From now on we shall focus on more interesting case 
of the repulsive Lorentz force (£ > 0). In this case there 
exist two qualitatively different types of the bounded mo- 
tion. Let p\ and pi be, as above, the minimal and the 
maximal values of p. If p2 < p* the right hand side of 
Eq. (|65[) is always positive. The motion has no curls, 
and 4> monotonically grows with time. This motion is 
modulated by oscillations in the radial direction. 

For pi > p* the motion of the particle is quite different . 
When the particle moves in the domain (pi, p*), we have 
d(j)/da > 0, and the angle <fi increases, whereas for the 
motion in the domain (p*,p2), we have d4>/da < 0, and 
the angle <f> decreases. The corresponding trajectory has 
curls. Increase of 0, corresponding to pG [pi, p„) , is not 
exactly compensated by its decrease, corresponding to 
p G {p*,P2\- As a result, there is a drift of the particle 
in the positive (^-direction. The critical type of motion 
corresponding to pi = p*, for which the trajectory is 
similar to a cycloid, is singled out by the condition £ = 
£*, where 



VI. BOUNDED TRAJECTORIES 

If the effective potential is a monotonically growing 
function of p, a charged particle always starts its mo- 
tion in the vicinity of the black hole, and after reaching 
the turning point, it falls down into the black hole. If 
the effective potential has extrema at p = p ma x,min, then 
for p < p max and £ 2 < U(p max ) a charged particle has 
similar motion. Here we shall not discuss such type of 
motion. Instead, we focus on study of bounded trajec- 
tories. For such trajectories the dimcnsionless radius p 
changes between its minimal p\ and maximal p- 2 values, 
Pmax < Pi < Pmm < P2, (see Figure EJ). These trajec- 
tories are more interesting for astrophysical applications, 
for example for analysis of the motion of particles in the 
black hole accretion disk. The specific energy £ for such 
trajectories obeys the relation 

^min ^ £ ^ ^max ■ (^^) 

The radial motion is periodic with the period 

dp 



Aoy = 2 



l>2 



, P1 ^£ 2 - U( P ) 
Let us consider equation (|4"Tj) for the angular variable 



(64) 



1/2 



(66) 



All the three types of bounded trajectories, with curls, 
the critical, and without curls, are schematically illus- 
trated in Figure O One can see that the bounded trajec- 
tories are similar to the trajectories of a charged particle 
moving in a weak gravitational field discussed in Sec- 
tion III. 

For a given b the motion can be specified by the con- 
served quantities £ and I. Different regions in the (£, /)- 
plane correspond to different types of motion. Let us 
discuss this in more detail. We assume that i > 0. The 
condition U tP = determines p ma x and p m in as functions 
of I. Substituting these functions into expression (|4"2")l 
for the effective potential one obtains the corresponding 
values = U(p maXtmin ,£). For the limiting value 



max,min 



I = 1+ , corresponding to the innermost stable circular or- 
bit, the two curves £ max (£) and £ m i n (£) meet each other. 
For I > £ + , the curve £ ma x(-0 is always above the curve 
£min(£)- At the points of the extrema of the function 
£ 2 (£) one has 
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FIG. 6. Types of a bounded trajectory corresponding to I > 0. Arrows illustrate direction of motion of a charged particle. 
Circular arcs represent the stable circular orbit defined by p = p m m + • (a): £+ > £,. Points A's correspond to p = pi, and point 
B corresponds to p = p2 > p*. Points C and D are turning points, where d<f>/da = 0. (b): £+ = £*. Points A's correspond to 
p — pi, and point B is a turning point, which corresponds to p — p2 = p», where d<f)/do — 0. (c): £+ G [£ m i n+ ,£,). Points A's 
correspond to p = pi , and point B corresponds to p = p2 < p* . 



where p = p 



min,max' 



Using Eq. (|52j) we conclude that 



d£ 



max, mm 



dt 



> 0. 



(68) 



Because corresponds to the positive slope of the po- 
tential U, one has £*(£) > £ m in(0- Figure [7] illustrates 
different types of motion for a given magnetic field B. 
Two curves on the 'energy-momentum' plane in this plot 
correspond to the maximum and the minimum of the 
effective potential. The curve £\ = £l represents the 
critical energy function of the angular momen- 

tum. The motion for the parameters below this line, in 
the domain /, is without curls, while the domain II cor- 
responds to trajectories with curls. 

Such a plot is convenient, for example for the discus- 
sion of the following problem: Consider a particle with 
the parameters E and L in the domain /. Let the parti- 
cle receives an additional portion of energy AE > and 
angular momentum AL. This process moves the point, 
representing the particle in {E, L) variables, to a new po- 
sition. If the point is moved to the domain II, then the 
excitation changes the trajectory of the particle, which 
becomes a curly one. 

Using expressions ([5]), (|59")l , and we can present 
the critical energy £ * in the dimensional form as follows: 



E 1 



2L 



> 0. 



(69) 



This expression establishes a relation between the param- 
eters E, L, and B, corresponding to the critical motion. 
For given B and L the motion is critical if E = . For 



the critical energy E* the bounded trajectory has cusps, 
one of which is illustrated by point B in Fig. [51(b). One 
can use Eq. to express the magnetic field B for the 
critical motion in terms of the energy E and the angular 
momentum L 



B,„ 



2L 

qr 2 g 



E v 



(70) 



If for fixed values of E and L the magnetic field is larger 
than B*, the motion is curly. Thus, keeping the energy 
and angular momentum of a charged particle fixed, one 
can change the type of its motion by changing the value 
of the magnetic field. 



VII. APPROXIMATE SOLUTION FOR 
BOUNDED MOTION 

To analyze properties of bounded trajectories near the 
minimum of the effective potential £^ in = U(p m - m ) one 
can expand U as follows 



U = £ min +LV (p- p D 



(71) 



where us is defined by Eq. (|T7| . and the dots denote 
the omitted higher order terms of the Taylor expansion. 
In the linear approximation (|7ip , equation (|40[) takes the 
following form: 



dp i 

— ] =£ - £ min - u) {p - p n 



(72) 
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hole in the counter-clockwise direction with the average 
angular velocity equal to Po = ^P m j n — b. This motion 
is modulated by the radial oscillations of the frequency 
w . Combination of the radial and azimuthal oscillations 
with the drift results in trajectories, qualitatively similar 
to those shown in Figs. [5] and [51 The critical solution 
is singled out by the condition that the trajectory has 
cusps. It happens for the following amplitude: 



A 



Po 
Pi 



2pt 



■{p 



2 2 

rmin 



(79) 



Since — p m - ln <C p m in, this condition implies that A w 
P* — Pmin or £ sa £*, which corresponds to the exactly 
calculated bounded trajectory with cusps. 

The average (drift) velocity of the guiding center in the 
direction of the increase of 6 is 



of b. Point A, where the curves fm ax (^) and f^ lin (^) join, 
corresponds to the innermost stable circular orbit. In the 
domain / between these curves bounded trajectories have no 
curls. In the domain // between the curves El and ^ ax 
the trajectories are of the curly-type. They degenerate into 
unstable circular orbit at point B. Bounded motion occurs in 
the regions / and II. This particular plot is constructed for 
b = 1/2. The qualitative behavior for different values of b is 
the same. 



V = PminA) 



Pn 



bpn 



(80) 



The same velocity defined with respect to the time t, 
differs from v by the Lorentz factor t given in (|20p . and 
is of the form 



V 



111111 



— 0*-/4m)- 



(81) 



For large values of b (see Eqs. © and (|7J|) equation 
U„ = can be solved approximately as follows: 



Analogous expansion of the angular velocity (|4Tj) near 
P = Pmin gives the linearized equation for <fi 



— = Po + Pl(P~ Pmin) , 



where 



Po 



b, Pi 



2t 



(73) 



(74) 



The validity of the linearized approximation requires that 



\P - Pmin | < 



|C^,ppp(Pmin)| 



, \p - Pmin | < Pmin ■ (75) 



We assume that these conditions are satisfied. 

Integrating Eqs. (|72[) and ([75} we derive an approxi- 
mate solution for a bounded trajectory of a charged par- 
ticle 



p(a) = p min + A cos(w cr) , 
P a H : 

Jf2 _ C2 
V ° °mh 



4>{u) = Po<J + ^— sin(cj a) , 
u>„ 



(76) 
(77) 

(78) 



corresponding to the initial conditions p(0) = pi = 
Pmin + A and 0(0) = 0. For £ > this solution de- 
scribes the motion of a charged particle around the black 



1 



lb»l 



86V* (P* - 1) 



(82) 



According to this expression, if £ is fixed, then p m i n is 
a decreasing function of b. Thus, strong magnetic field 
shifts stable circular orbits toward the black hole horizon. 
For b 3> 1 in the approximation of a weak gravitational 



field, that is when p n 
velocity is 

1 

V 



> 1 and E 



Abp\ 



min 



2qBr 2 



9o_ 



the drift 



(83) 



where r Q = p mm r g and go = r g /(2rg). This expression 
coincides (as expected) with expression ([36)) for £ w 1. 

For b S> 1 in the approximation of a strong gravita- 
tional field, that is p m i n ~ 1, one has 



V oc 



' Sir. 



(84) 



where fl 2 ^ is the Keplerian frequency (U]) calculated for 



Pn 



The proportionality constant in ([843 



depends on the value of the angular momentum L and 
is of the order of unity. We see that in both the cases 
of strong and weak gravitational fields the gravitational 
drift velocity is proportional to B^ 1 . 

Another quantity which characterizes a bounded tra- 
jectory corresponding to the repulsive Lorentz force is 
the ratio N of the frequency of the radial oscillations uiq 
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to the average angular velocity (3 Q . Using Eqs. (|48|) and 
(f?4")) we obtain 

N _ ^0 = V >i in (3pmin ~ 1) + P^Pmin ~ 1) 

~~ fl> /9 1/2 (rt 2 - P 2 ■ ) 

rminvr* rmin/ 

For A > — /3o / /3i this ratio gives the number of curls per 
one revolution of a particle around the black hole. This 
ratio depends on p m ; n and p*. Using the approximation 
((82)1 for b > 1 we have 

N\ b>>1 ^8b 2 pl /2 (p*-l)V 2 . (86) 

Thus, the stronger magnetic field, the greater is the ratio 
and the number of curls. 



VIII. SUMMARY 

In this paper we studied motion of charged parti- 
cles in the equatorial plane of a weakly magnetized 
Schwarzschild black hole. We analyzed properties of the 
corresponding effective potential due to the combined 
gravitational and Lorentz forces acting on a charged par- 
ticle. We gave a simple analytical proof that this poten- 
tial either has two extremal points or none in the black 
hole exterior. The critical case, when these extrema co- 
incide, determines a position of the innermost stable cir- 
cular orbit (see also 0-[ll|). We obtained expressions 
for the radii and the corresponding angular momenta of 
the innermost stable orbits in the approximation of the 



strong magnetic field (6 > 1). A similar expression for 
the radii in the limit b — > oo was given before in [ill [l3| . 

In our analysis we mainly focused on the study of 
bounded trajectories of charged particles. Such trajec- 
tories may be considered as an approximation to charged 
particles motion in the black hole accretion disk, when 
their mutual interaction is neglected. We constructed an 
approximate solution to the dynamical equations which 
represents a bounded trajectory localized near the stable 
circular orbit. 

As a result of our study, we found that if the Lorentz 
force acting on a charged particle is repulsive, its bounded 
trajectory can be of two different types: with curls and 
without them. The critical trajectory, which has cusps, 
separates these two cases. We calculated the critical 
value of the magnetic field for the critical trajectory. We 
calculated also the number of curls per one revolution and 
found that its maximal value grows with the increase of 
the magnetic field. Using the approximate solution we 
found the gravitational drift velocity of the guiding cen- 
ter of a bounded particle trajectory. 
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